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Abstract. It is assumed that Be star discs are driven by 
viscosity. Emission from disc models is calculated and is 
confronted with continuum observations. It is found that 
the outflowing viscous disc models can reproduce the ob- 
served IR continuum emission. However, to exist as out- 
flowing discs, either the discs are significantly acted upon 
by the stellar radiation field and/or there is significant 
cooling with radius in the disc. The energy generated via 
viscous dissipation is calculated and shown to play only a 
minor role in the energy balance of the disc. 

A scenario whereby a B star may change into a Be 
star (and vice versa) by generating (reaccreting) the disc 
is suggested. 

Key words: stars: emission-line. Be: mass loss: rotation 
- circumstellar matter 



1. Introduction 

Be stars are fast rotating stars (Slettebak 1982) with rota- 
tion rates up to break-up, and a modal value of ~ 70% of 
break-up rotation (Porter 1996). The circumstellar mat- 
ter around Be stars is thought to consist of two distinct 
regions: a diffuse polar stellar wind and a dense equato- 
rial "disc" (Dachs 1987, Slettebak 1988). The nature of 
the polar wind seems to be well understood in the con- 
text of radiation-driven wind theory (Castor, Abbott & 
Klein 1975, Kudritzki et al. 1989). However, the mecha- 
nism generating the equatorial disc is more elusive: several 
theories have been proposed including wind compression 
(Bjorkman & Cassinelli 1993), outflowing viscous discs 
(Lee et al. 1991), and wind bistability (Lamers & Paul- 
drach 1991). Recently it has been found that stars with 
high rotation will generate a metallicity enhancement in 
their equatorial planes (Porter 1999), which through the 
metallicity dependence of the radiation driving parame- 
ters may change the nature of the wind there. 
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Wind compression has already been shown to be insuf- 
ficient on theoretical (Owocki, Cranmer & Gayley 1996) 
and observational (Porter 1997) grounds. The major prob- 
lem with wind bistability is that the disc appears to be 
rotating in a Keplerian fashion (e.g. Dachs et al. 1986, 
Hanuschik 1989, 1996), whereas any wind causing the 
equatorial enhancement is likely to conserve angular mo- 
mentum, and should not be observed to be rotating faster 
than the star - Hanuschik (1996) has shown that the half- 
line widths of emission lines created in the disc are indeed 
larger than v sini of the star. 

The theory which seems to be most applicable to Be 
star discs is that of viscously driven outflow (Lee et al. 
1991). Here angular momentum is added to the inner edge 
of the disc increasing its angular velocity to slightly super- 
Keplerian. The disc interacts with itself via viscosity and 
angular momentum is transported outwards. The angular 
momentum source has been suggested by Osaki (1986) to 
be due to non-radial pulsations dissipating in the atmo- 
sphere of the star. The outflowing viscous disc model has 
several circumstantial corroborations - the variation in the 
line asymmetry (the V/R variation) can be naturally gen- 
erated in Keplerian discs (Papaloizou et al. 1992, Okazaki 
1991). Excess IR emission generated by the disc is seen 
to disappear over a timescale of months (Hanuschik et al. 
1993) which is similar to the viscous timescale in these 
discs ~ this is perhaps indicative of phases of outflow and 
inflow within the disc. 

Although the outflowing viscous disc model appears to 
be the most likely, it has never been subjected to observa- 
tional tests - can these discs generate the IR excess? What 
is the source of viscosity? How and why do they dissipate 
and regenerate? In this investigation. Be star discs are as- 
sumed to be outflowing viscous discs and an attempt is 
made to answer these questions. 

In §2 aspects of the outflowing viscous disc model are 
considered along with their application to Be star discs. 
In §3 the energy balance within the disc is studied, and 
the theory is applied to a test case in §4. Discussions and 
conclusions are presented in §5 and §6. 
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2. Can outflowing viscous discs exist? 

2.1. Viscous discs 

The gas dynamics of viscous discs has been analysed and 
developed by several authors (Shakura & Sunyaev 1973, 
Pringle 1981, Frank, King & Raine 1992). Here, some com- 
ments are made regarding the structure of such discs (fol- 
lowing arguments of the references above). First, a disc 
interacting with itself via viscous stresses transports an- 
gular momentum outward. In order for viscous stresses to 
exist in a disc, there must be (i) a non-zero viscosity and 
(ii) some differential rotation in the disc. Conservation 
of angular momentum in a cylindrical co-ordinate system 
iR,(t>,z), yields 



4(^«^"> + l^(«^-«'">^^i 



(1) 



where S is the surface density of the disc, fl is the an- 
gular velocity, and vr is the radial "drift" velocity. The 
viscous torque is G = 2TTRvYjR^{dQ./dR) where v is the 
viscosity. If a ring of matter is placed in a Keplerian orbit 
and allowed to interact with itself through viscous stresses 
then part of it will move to smaller radii having w^ < 
and parts of it will move outward w/j > (e.g. Pringle's 
1981, fig.l). This is solely due to the redistribution of an- 
gular momentum in the gas. The outflowing (inflowing) 
parts have increased (decreased) angular momentum with 
respect to the initial state. Time steady {d/dt = 0) solu- 
tions which show no drift velocities w^ = occur when 
G is a constant with radius, i.e. dVt/dR oc {vYiR^)~^ , al- 
though these solutions have limited physical application 
in discs around stars. These static discs still transport an- 
gular momentum outwards, even though the gas in them 
is not moving radially. 

For an accretion disc, matter is added at the outer part 
of the disc. The viscous stresses lower the angular momen- 
tum of a ring of gas which in response moves inward and 
increases its angular velocity until it it regains rotational 
support. This mechanism creates the accretion flow: the 
gas moves inward vr < ^ throughout the disc. The disc's 
velocity field is typically rotationally dominated with ro- 
tational velocities very close to Keplerian, and a subsonic 
radial velocity. 

Eq.l also permits solutions with positive radial drift 
velocity vr everywhere i.e. outflowing discs. This sort of 
disc requires an angular momentum source at the inner 
boundary - the disc-star interface - as well as a gas source 
taken to be the atmosphere. The disc recieves angular mo- 
mentum from the star which will itself spin down. It is 
exactly this sort of model which has been suggested by 
Lee et al. (1991) to apply to the discs of Be star and are 
investigated below. Porter (1998a) has already condisered 
the angular momentum evolution of Be stars for this sort 
of disc, and finds that the observational result that there 
is little or no angular momentum evolution of Be stars 
(Zorec & Briot 1997, Steele 1999) may be explained if the 



viscous stresses at the inner edge of the disc are small or 
that the disc is present intermittently in the lifetime of the 
star. A small couple applied from the star at the inner edge 
leads to slow outfiow velocities of vr <^ 0.01 km s^^at the 
inner edge of the disc. Note that an outflowing disc cannot 
exist if there is no couple from the star at the inner edge. 
Aside from Lee et al.^s (1991) initial suggestion and 
modelling of Be star discs, outflowing viscous discs have 
been discussed by Pringle (1991), Narita, Kiguchi & 
Hayashi (1994) and Okazaki (1997). Narita et al.^s numer- 
ical analysis showed that if the central object rotates fast 
enough, then the disc docs indeed become an outflowing 
one. Let us assume then that Be star discs are outflow- 
ing viscous discs and subject them to some observational 
tests. 

2.2. Angular momentum transport in a viscous disc 

The existence of an outfiowing disc is dependent on angu- 
lar momentum being added at the inner edge of the disc 
whereas the dynamics of the disc is determined by the vis- 
cous stresses: if only this initial couple is present then a 
ring of gas will move away from the star until it finds its 
Keplerian orbit, at which point it will stop. The evolution 
of the disc away from the inner boundary is governed by 
the viscosity. 

The rotational velocity of a outflowing viscous disc is 
v^ « y^GM^/R, where M^, is the mass of the star, and 
G is the gravitational constant. Assuming that the disc is 
steady (so d/dt = 0), the equation for transport of angular 
momentum becomes 
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where S is the surface density of the disc, vr is the radial 
velocity, and v is the viscosity. 

The surface density is S = p{R,0)H, where p{R,0) 
is the density in the equatorial plane at radius R, and 
H — R Cs / Vcf, is the density scale height (c^ is the sound 
speed). The disc is isothermal in the z direction and as- 
suming that the density in the equatorial plane p{R, 0) is 
a power law, the density field is 



p{R, z) ^ po[—- 



R 

R. 



exp 



2iJ2 



(3) 



where po is the density at the inner boundary. Note that 
the numerical simulations of outflowing discs by Narita et 
al. (1994) find that steady state outflowing discs do indeed 
have power law surface densities from R > 2i?, , justifying 
the power-law assumption in eq.3. 

Following convention, an alpha prescription is used for 
the viscosity ly = aCgH (Shakura & Sunyaev, 1973). Fi- 
nally, to ensure the disc is as general as possible, the tem- 
perature of the disc is assumed to follow a power law: 
Td = O.SToff (i?/i?, )"", where T^s is the effective tem- 
perature of the star. 
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2.3. Inflow or outflow? 

Let us assume that we have a disc which is viscously in- 
teracting with itself. Is it possible to tell whether the gas 
flow is inward, outward or zero? In asking this question, 
the nature of the addition of torque to the inner boundary 
of the disc for outflowing discs, or that of the angular mo- 
mentum removal for inflowing discs has been ignored. An 
attempt is made to ascertain whether outflowing viscous 
discs are credible Be star disc candidates. 

For vr > then the radial exponent of the term in 
brackets in eq.2 must be less than zero, ensuring the differ- 
ential is negative and hence the right hand side is positive. 
Inserting eq.3 and the a viscosity into eq.2 and collecting 
powers of radius R together, then the right hand side is 
proportional to —d / ^R{R^^^^"■^^^^'^) . Therefore the disc 
is an outflowing disc if 2n + 3to > 7, assuming that a is 
a constant (see later). It is worthy of note that when out- 
flowing viscous discs have been modelled, this limit on n 
and m has always been obeyed (e.g. see Narita, Kiguchi 
& Hayashi 1994 and Okazaki 1997). 

2.4. Confrontation with observations 

There does appear to be a problem with the outflowing 
viscous disc paradigm straight away when confronting it 
with observations. Cote & Waters (1987) found that the 
IR emission (originating in the inner 10s of i?, of the disc) 
is well fit with an isothermal disc (to — 0) having a den- 
sity power law index of 2 < n < 3.5. Values of n are larger 
when fitting the emission in the radio regime, correspond- 
ing to larger radii of the disc (Dougherty et at, 1991). 
According to the above work, the inner parts of the discs 
cannot be outflowing viscous discs, and must actually be 
accretion discs! 

If the outflowing viscous disc model is applicable to Be 
star discs then this discrepancy must be examined. First, 
consider the emission fitting procedure: typically both a 
constant opening angle for the disc and an isothermal disc 
are used. The constant opening angle leads to H (x R a.s 
opposed to H oc R^/'^ for viscous discs above. Inserting 
this slight modification into the angular momentum con- 
servation expression, leads to the new limit 2n + 3m > 6 
for outflow. 

However, comparing this with the derived values of n 
does not really change the conclusion that almost all discs 
should be accreting! We must be wary of taking the fits 
at face value - usually the disc was assumed a priori to 
be isothermal, and the parameter m was not allowed to 
vary. If a non-isothermal disc is used then the fits to the 
IR excess may change. 

What can be gleaned from this apparent paradox that 
all Be star discs should be inflowing, given that there is a 
strong belief that they are in fact outflowing? It is unlikely 
that the fitting procedure produces values of n consistently 
too small in almost every cases (see §4). The conclusion. 



therefore, is that the outflowing viscous disc can not fit 
the observations, and cannot be applicable to Be stars. 

This is too hasty: aspects of the problem have yet to be 
considered. The first is the radiation field from the star: 
the total electron-scattering optical depth in the equato- 
rial plane is 



R. 
0.24 
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(JePoR* 



(4) 
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where ae — 0.35cm^g^^ is the electron scattering cross- 
section. Inserting typical values of n = 2.5, po = lO^^^g 
cm"'^ and R^ — 5Rq yields an optical depth of r = 0.5. 
Therefore the disc is optically thin to electron scattering, 
and hence the disc may be acted upon by optically thin 
lines, as suggested by Chen & Marlborough (1994) (note 
that the optical depth argument derived in §4.1 of Chen & 
Marlborough where they rule out radially subsonic discs 
is misleading due to their choice of wind parameters: vis- 
cous discs may have mass-loss rates of IO^^^Mq yr^^ with 
small radial velocities and still be dense enough to produce 
the IR excess emission e.g. Okazaki 1997). 

If regions of the disc are partially supported and driven 
outwards by radiation, then they will deviate from Ke- 
plerian rotation. This is because an insufficient amount 
of angular momentum is given to a ring of material as it 
moves outward. Hence the rotation velocity of the disc will 
lie between Keplerian (w^ oc r~^'^) and angular momen- 
tum conserving {v^f, oc r~^). This difference may be very 
difficult to discriminate observationally, and so radiation- 
driven discs cannot be ruled out. However, it should also 
be noted that the inclusion of optically thin radiation driv- 
ing does not necessarily produce non-Keplcrian discs (e.g. 
see Okazaki 1997). 

It is noted that the viscosity parameter a may be de- 
pendent on radius. The parameter a represents the dif- 
ference between the actual value of (size x velocity) of the 
turbulent eddies and the characteristic product of scale 
and speed in the disc (disc height x sound speed). In fact 
if the turbulence becomes supersonic then it is possible 
that a > 1 . In the situation currently considered, a dense 
rotationally-dominated disc is adjacent to the fast wind 
which is dominated by the radial component of velocity. 
In the interaction region there are large shearing velocities 
and therefore the interaction region will become unstable 
to Kelvin-Helmholtz instabilities. 

Consequently it might be expected that the typical 
turbulent velocity will be a function of the Mach number 
of the fast wind. In fact, if the eddy size is the typical scale 
height of the disc, then the maximum value that is possible 
if this is the case will be the fast wind Mach number - the 
wind typically reaches velocities of >1000 km s"-'^ after a 
few R^ , and with the typical temperature of the disc of 
10 km s^^, yield Mach numbers of A^ ^ 100. Therefore, it 
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is conceviable that may become larger than unity at quite 
a small radius, although a < A4 ^ 100. 

As the fast wind's velocity increases with distance from 
the star, then a might be expected to increase with radius. 
This aspect of the viscosity is difficult to examine and a 
full study is left to a future paper. It is expected that if this 
Kelvin-Helmholtz process is significant in determining a, 
then it should increase with radius. This, however, makes 
the gas in the observed discs more likely to be inflowing 
rather than outflowing. 

It is also noted that the disc may be outflowing if there 
is a sufficient radial temperature gradient in the disc (cor- 
responding to large to) . This requires the disc to cool sig- 
nificantly as it drifts away from the star. However, recent 
work by Millar & Marlborough (1998) casts doubt on a 
such a large temperature gradient being present. 



2.5. Accretion discs? 

The above paragraphs attempt to show that outflowing 
viscous discs are acceptable models for Be star discs. How- 
ever, to make that discussion the primary one for the rest 
of this paper, we must answer the question "why are Be 
star discs not normal accretion discs?" 

There are two prerequisites for them to be accretion 
discs: there must be a supply of gas to make the disc, and 
the gas must have larger specific angular momentum than 
Keplerian at the disc surface. The first of these two as- 
pects can be fulfilled at least for some stars with weak 
radiatively driven winds: Porter & Skouza (1999) have 
examined models of radiatively driven winds where the 
gas and radiation field decouple due to ion stripping (e.g. 
Springmann & Pauldrach 1992), before the flow becomes 
unbound from the star. The gas then stalls in the star's 
gravitational potential and reaccretes. This may provide 
the material to make an accretion disc if the shell then was 
concentrated in the equatorial plane, but it will not have 
enough specific angular momentum to form an accretion 
disc. If the gas conserves angular momentum (almost cer- 
tainly true for the outflowing phase), then when the gas 
reaccretes, it will have the same angular velocity as it left 
the star, i.e. sub-Keplerian. Consequently, unless there is 
some way of adding angular momentum to the wind as it 
stalls, this does not seem a feasible way to make accretion 
discs. A second possible reservoir of gas is a companion 
star, which is losing mass to the Be star - this does not 
have the angular momentum problem, although it does 
mean that all Be stars are in binary systems. 

On the balance of these points, it seems unlikely that 
Be stars discs are accretion discs because (i) not all Be 
stars are in interacting binary systems with the secondary 
star as the mass donor, and (ii) an accretion disc will not 
form when a star is reaccreting its own wind. 



3. Energy balance within the disc 

3.1. Liberation of energy via viscous dissipation 

The viscous stresses enables the disc to flow outwards and 
in doing so liberates energy. What happens to this energy? 
The timescale over which the viscous disc can evolve sig- 
nificantly Ty is defined: 
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where T4 is the disc temperature in lO^K. Inserting typical 
parameters leads to r^, '^ lO^s of days. Given that Be star 
discs are observed to be present for many years (a pre- 
requisite for them to be able to develop and evolve V/R 
variations in the lines), then a typical disc is stable over 
a viscous timescale. Therefore any energy liberated via 
viscous stresses must emerge as radiation. 
The luminosity emitted is 



dL 
dR 



97r ^ f R 
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(6) 



(Frank, King & Raine, 1992, p73), where the expressions 
for H , S and v have been inserted to obtain the second 
equality. This luminosity source is termed the "viscous 
luminosity" as it is the energy emitted caused solely by the 
viscous processes leading to outflow. This expression only 
gives the luminosity as a function of radius, and contains 
no spectral information. As the disc is optically thin in 
the z direction, then the standard approach to the disc 
emission (e.g. Section 5.5 of Frank, King & Raine 1992) 
are inapplicable, and an alternative discussion is needed. 

3.2. IR free- free and free-bound emission 

The continuum emission of Be star discs is limited to the 
IR-radio spectral regions (e.g. see fig. 11 of Poeckert & 
Marlborough 1978) and was identified as free- free and free- 
bound emission from an ionized plasma by Gehrz et al. 
(1974). The most successful empirical model used to cal- 
culate the emission is due to Waters (1986), and represents 
the disc using a density power law, and an opening angle. 
Three parameters determine the emission: the opening an- 
gle, the density at the star-disc boundary, and the density 
power-law exponent. 

The above expressions for a outflowing viscous disc 
changes the calculation of disc emission and so some of 
Waters' expressions are now restated. The optical depth 
T^ of the disc may be calculated at a frequency v. The 
expression is simpler than Waters' expression as the inte- 
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gral through the disc (equivalent to his C{9,rL)) may be 
completed analytically. Using Waters' notation 



tAR) = Xx X,d i?-2»+3/2-™/2 
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where A is the wavelength in cm, gff and gfi, are the Gaunt 
factors for free-free and free-bound emission respectively, 
z"^ is the mean squared atomic charge, 7 is the ratio of the 
number of electrons to ions, /i is the mean atomic mass, 
and Csfi is the sound speed at the star. 

The intensity Ii, (R) is calculated in the same way as in 
Waters (1986), and hence the spectrum of the disc may be 
calculated. However, here it is the total emission at a given 
radius which provides a direct comparison with outflowing 
viscous disc models. The differential energy liberated is 
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(equivalent to Waters' eq.l2). Finally to obtain an expres- 
sion to relate to eq.6, this is integrated over all frequencies: 



dL 9 

dR - ^'^^ 



h{R)dv. 



(9) 



3.3. How important is the viscous luminosity? 



There are now two expressions for differential disc lumi- 
nosity deriving from the viscous stresses (eq.6) and the 
free- free, free-bound emission of the plasma (eq.9). Here 
the fraction of viscous luminosity to the free-free, free- 
bound emission of the disc is considered. This is an im- 
portant check on the model - if the viscous luminosity 
(eq.6) is larger then the disc emission (eq.9) then the disc 
cannot be driven by viscosity. Eq.8 provides a prescription 
for the disc emission which may be directly related to the 
observations by integrating over radius to produce a spec- 
trum, whereas eq.6 provides the theoretical emission due 
to viscous dissipation. 

The density parameter po and the power law exponent 
n, can be constrained by fitting the spectrum of the disc 
with the spectrum produced by Waters' method. The ra- 
tio T of luminosity emitted from free-free and free-bound 
transitions to the viscously dissipated energy is 



(10) 



The most poorly constrained parameter is then the vis- 
cosity parameter a which is a priori limited to the region 



Table 1. Best fit parameters for the IR excess of x-Oph. Model 



1 is isothermal and so m = 0. 






log(po) (g cm"^) n 


m 


2n + 3m 


Model 1 -11.22 2.20 


(0) 


4.40 


Model 2 -11.24 1.90 


0.23 


4.49 



< a < 1 . However this may an underestimate if the tur- 
bulence becomes supersonic (see §2.4). For the current cal- 
culation a is set to unity, although the limit that a ^ 100 
is kept in mind. 



4. A case study x-Oph 

To illustrate the discussion above, a test case is considered 
- this provides us with a "real" example. As an example, 
the observations of the B2IV star x-Oph are used (Waters 
1986). The star was assumed to have a mass, radius and 
temperature of lO.OM© , 5.7i?, , and 22,500K. 



J^.l. IR continuum excess fits 

To constrain the power law exponents, the IR excess is 
first fitted by the outflowing viscous disc outlined in §2 
using the procedure due to Waters (1986). Two fits to the 
excess were calculated: one assumed an isothermal disc 
(Model 1) and the other allowed the temperature power 
law exponent m to vary (Model 2). The disc temperature 
at the inner edge was assumed to be 18000K (i.e. 0.8 of 
the effective temperature). In this latter case, the emis- 
sion was calculated only for radii which had temperatures 
larger than 10'*K to ensure that the region was ionized. 
In both cases the disc is assumed to have a large radius 
(i?disc > 50i?* ). The best fit parameters for the models 
are presented in table 1, and are shown in fig.l. 

Note that the derived parameters are not the same as 
those of Waters (1986) as the disc model used in that work 
differs from that used here. Model 2 is formally a better 
fit to the data - the reduced x^ values are 1.1 and 0.7 for 
Models 1 and 2 respectively. Although Model 2 fits the 
data better, it's reduced x^ value of 0.7 indicates that it 
is possibly "over fit" and so it does not provide strong ev- 
idence for non-isothermality in the disc. The temperature 
falls below lO^K at i? ss 13i?* providing a natural outer 
boundary to the emitting part of the disc ~ this radius 
is low when compared to fits including millimeter obser- 
vations of other stars (Waters et al. 1991). Millimeter- 
centimeter data provide constraints on the extent of the 
outer parts of the disc, and hence the temperature expo- 
nent TO. Column 5 of table 1 shows the combination of 
2n -|- Sto, which must be larger than 7 for a non-radiation 
driven outflowing viscous disc - both models fall signifi- 
cantly short. 
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Fig. 1. The IR continuum excess emission as a function of 
wavelength. Z — 1 is the fractional emission excess over the 
stellar photosphere. All data is taken from Waters (1986). The 
solid line is the fit from the isothermal Model 1 , and the dotted 
line corresponds to the non-isothermal Model 2. 



4-2. Viscous luminosity versus free-free & free-bound 
emission 

Now the power law exponents have been fitted from the 
spectrum, the relative contribution to the total emission 
of the viscously dissipated energy can be calculated. The 
ratio J- is evaluated from eq.lO and is shown in fig. 2 for 
both models. 

Clearly the free-free and free-bound emission domi- 
nates over the viscously dissipated energy, by of order 1000 
for both models. This result states that the energy balance 
within the disc is not dominated by the energy produced 
by the shear motions within the disc itself, and hence must 
be dominated by the stellar radiation field. Even with the 
maximum value of a '-~^ 100, the free-free and free-bound 
emission still dominates, making the result secure. 



5. Discussion 

The previous sections have provided a consistent descrip- 
tion of Be star discs as outfiowing viscous discs. This study 
links outflowing viscous disc theory with observations and 
finds that the two are indeed compatible (the first out- 
flowing viscous disc paper by Lee et al. 1991 produces 
very high disc densities due to their choice of high disc 
mass-loss rate which can be ruled out by examining the 



Fig. 2. Ratio of free-free and free-bound disc emission to 
the viscous emission from eq.lO. The solid line corresponds 
to Model 1, the dotted line to Model 2. 



IR continuum emission). This study therefore adds to the 
growing amount of evidence that Be star discs are in re- 
ality viscously driven. There is already substantial work 
regarding instabilities in Keplerian discs which give rise 
to density perturbations which produce asymmetric line 
profiles (e.g. Okazaki 1991). Viscous discs can now explain 
almost all of the observations of Be star discs. 

One major observational aspect of Be stars is that they 
change phase from Be-Bc-shell-normal B star, occurring 
apparently at random (although the disc may disperse and 
reappear in a more orderly fashion e.g. /iCen, Hanuschik 
et al. 1993). How may the outflowing viscous disc model 
explain this? There may be a clue to these changes in the 
suggestions made in §2.3 regarding the regimes where the 
observed discs can be outflowing discs. 

Two possibilities have been speculated upon: (i) if the 
disc is supported by radiation, then it should be noted 
that the electron scattering optical depth is typically in 
the range 0.1-1.0 (eq.5). With this so close to unity, vari- 
ations in the density at the star-disc boundary of factors 
of a few may be enough to decrease the radiative driv- 
ing through the disc. This loss of radiative support will 
lead naturally to an accretion phase where the disc falls 
back on to the star. The density variations at the star-disc 
boundary could, in principle be created as a by product 
of disc warping (Porter 1998b), or non-radial pulsations of 
the star (e.g. Lee & Saio 1990). Alternatively, (ii) the disc 
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may be a outflowing viscous disc due to a possible large ra- 
dial temperature gradient. However the recent modelling 
of Millar & Marlborough appears to rule this out. 

With the outflowing viscous disc model able to explain 
the current observations it is pertinent to ask whether fu- 
ture observations could confirm the paradigm. The cur- 
rent major weakness with the theory regards the input of 
angular momentum at the inner boundary. Although this 
process has been suggested to be dissipation of non-radial 
pulsations in the atmosphere (Osaki 1986), it has not been 
clearly demonstrated. If this is the underlying mechanism, 
then there should be statistical correlations of disc and 
pulsation parameters - an observational aspect which is 
still to be resolved convincingly. 

6. Conclusion 

This paper has produced several points: first, outflowing 
viscous discs can only exist around Be stars if the disc is 
partly driven outwards by the stellar radiation field This 
has been discussed and shown to be viable for actual Be 
star discs. 

Secondly, it is found (confirming a priori expectations) 
that a viscous disc can account for the observed excess IR 
emission of Be stars. Also, it has been found that the en- 
ergy balance in the disc is dominated by the stellar radia- 
tion field. The energy liberated in viscous dissipation has 
been calculated and is shown to be a small fraction of the 
observed luminosity of the disc (^ 10~'^-10~^). These sec- 
ond two points are the first time it had been demonstrated 
that outflowing viscous discs are energetically allowable 
models. 

The underlying reason why Be stars go through phases 
where their discs are lost has been speculated upon. This 
may be due to blocking of radiative support and driving of 
the disc at the inner regions - causing the disc to change 
to an accretion disc. 
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